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Various types of mixed spin two-dimensional Heisenberg networks are investigated by means 
of Monte Carlo simulations. This study aims at interpreting quantitatively the thermodynamical 
properties of two-dimensional molecule-based magnets recently synthesized. The proposed model 
requires that: (i) one of the two magnetic centers has a spin large enough to be treated as a 
classical spin; (ii) the zero field Hamiltonian is isotropic; (iii) the quantum spins have only classical 
' spins as neighbours. The quantum Hamiltonian is then replaced by a classical one with effective 

ferromagnetic interactions. The temperature dependence of both the specific heat and magnetic 
susceptibility are calculated. The effect of the lattice geometry is analysed. 
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I. INTRODUCTION 



CO 

A rather large number of molecule-based magnets have been synthesized and investigated in the last few years 

ajl],D. They correspond to low-dimensional magnetic systems, either quasi-one-dimensional or more recently, 

quasi-two-dimensional ]9|Jl0[|. The one-dimensional compounds are well modeled for equally-spaced magnetic chains, 
alternating chains involving a unique kind of spin carrier and two kinds of interaction pathways, mixed spin chains, 
£h ' ladder- type double chains, etc... 

This paper is devoted to two-dimensional Heisenberg mixed spin compounds in which one of the spins is large 
i^i - enough to be treated as a classical spin and the other is normally treated as a quantum spin. This work is motivated 
by the synthesis of novel two-dimensional magnetic materials. So far, two types of two-dimensional magnetic lattices 
t-H , have been described. Both have a honeycomb-like structure. The former type is obtained using oxamate as bridging 
ligand. Mn 2+ ions in octahedral surroundings are located at the corners of the hexagons, and Cu 2+ ions in elongated 
tetragonal surroundings are located at the middle of the edges. A strong antiferromagnetic interaction is propagated 
between the Mn 2+ and Cu 2+ ions through the oxamate bridge, so that the intralayer interaction is very large as 
compared to the interlayer interaction. In these compounds the layers are negatively charged, and the nature and 
• magnitude of the interlayer interaction is governed by the size of the countercations situated between the layers. In 
f-*. | the case of the NBu 4+ countercation a long range magnetic transition was observed at 15 K, probably due to the 
synergy between a very weak magnetic anisotropy and a ferromagnetic interlayer interaction |11[] . It has been shown 
in a previous paper |T^ ] that our model leads to an excellent interpretation of the magnetic properties of this material. 
The parameters were found as J — 47.6 K, gun = 2.0 and gc u = 2.2. These values are very close to those obtained for 
both Cu 2+ Mn 2+ pairs || and chains involving the same bridge. The value of the interaction parameter essentially 
depends on the nature of the bridging network, and is not very sensitive to the spin geometry. 

The latter type of lattice is realised in a series of two-dimensional oxalate-bridged bimetallic compounds which have 
just been synthesized [|l0]]. In that case the two kinds of magnetic ions alternate at the corners of the hexagons. The 
O ■ layers which are again negatively charged are separated by countercations. The general formula of these componds 
is'(NBu 4 )[M // Ru /// (ox) 3 ] where M /7 is either Mn" (S =5/2) or Fe" (S =2) or Cu 11 (S =1/2) and the spin carried 
by the high-field Ru 3+ ion is 5(Ru) = 1/2. The interaction is ferromagnetic for Mn and antiferromagnetic for Cu 
and Fe. In the latter case a ferromagnetic transition occurs at T c = 13K. Let us mention that other compounds with 
the same structure have been reported |p3|d!6|. The nature of the spin carriers, however, does not allow to use the 
■ - - ■ classical-quantum spin approach. 
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LATTICE (a) LATTICE (c) 




FIG. 1. For each cases, on the left hand side, the original lattice, where the squares stand for quantum spins, and the circles 
for classical spins. On the right hand side, the effective lattice of classical spins, with in bold the unit cell. 

The aim of this paper is to study the influence of the geometry of the spin lattice on the thermodynamical properties 
of such systems. In this way, we derive from the mixed quantum-classical Heisenberg model a purely classical one 
with an effective interaction which depends on the lattice geometry. This model is analysed by means of Monte Carlo 
(MC) simulations. We investigate the different realisations of the hexagonal lattice described above, but also two 
configurations based on the square lattice. 

The paper is organized as follows: the model is developed in the first section, then the Monte Carlo analysis is 
presented, and the results are discussed in the last section. 

II. THE MODEL 

Let us write down the Heisenberg spin Hamiltonian as: 

N c N Q 

{ij) j=i 1=1 

Here Sj is the large spin operator (5/2 for Mn, 2 for Fe) which will be approximated by a classical vector, Ss 
where s is a unit vector and S = y/S^(S^ + l).The 1/2 spin quantum operator (for Cu or Ru) is denoted by 
— g<Tj with a - , the Pauli matrices. The interaction parameter J is positive for an antiferromagnctic interaction 
(in the following, we only consider this case); H is a weak magnetic field applied along the z direction. < ij > stands 
for a pair of nearest neighbour spins, Nc is the number of classical spins and Nq is the number of quantum spins. 

Two kinds of magnetic lattices are investigated, the hexagonal one, which has been realised experimentally, and 
the square one. For each lattice, the spins can be arranged in two fashions: either the classical spins are attached at 
each vertex, and the quantum spins occupy the middle of the links, or the classical and quantum spins alternate at 
the vertices of the lattice. The structures are schematized in figure 1 (a) to (d). 

In all cases, a quantum spin is surrounded only by classical ones. Therefore, the partition function can be factorized 
with respect to the quantum spin operators: 



Z(T,H) = 




In this expression, V(i) is the set of labels of the classical spins nearest neighbours of the quantum spin at site i. 
Let us call this set of classical spins a unit cell. Depending on the lattice, the unit cell is a link (lattices (a) and 
(c) of figure 1), a triangular plaquette (lattice (d)) or a square plaquette (lattice (b)). Due to the factorized form of 
Z(T,H), the quantum spin dependence can be traced out to give a fully classical partition function: 
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Z(T,H) = J \Jl dn A |n 2cosh 11 ~\ PJS g, 8 ' + \P92lx B He z || j exp (pg lf x B S H-J^. 



(1) 



where {r} is the set of unit cells on the lattice. 

The original mixed spin system is then equivalent to a classical one with an effective ferromagnetic interaction 
between the classical spins on the plaquettes {T} which becomes (in zero field): 



H e & = ~k B T V In ( 2 cosh || K £ Sj || ) 

{r} V ^ T ) 



In the following, we shall use the notations 
1 JS 



l k B l jgr 



W Z (T) = W(r) • e 2 



The various observables, like the heat capacity Cy = k B (3 2 -g^hiZ{T,Q) and the zero field suceptibility 



_ k B T d 2 i 7 



H=0 



are simple generalisations of the ones defined in ref. [|12| and can be expressed as ensemble 

averages with respect to the Boltzmann weight e _/3W<!ff /Z(T, 0). By defining E = -r\JSJ2{r} w ( t ) tanh(K W (T)) , 
we find that the internal energy and the specific heat are given by 



U = 
Cv = k B (3 2 



<^V«-< B &J with * = E 

{r} 



KW(T) 



n 2 



cosh(iT W(T)) 



The molar magnetic susceptibility is obtained from eq. ([!]): 

where iVjvf is the number of molecules and P, Q, and R have the following expression^ 



(2) 



(3) 



\i=l / 



(S 81 ) * fe> ,r >) 



a= E s "*( r ) -E^^+El^-^ 11 )) 



jr} / {r} {r} 
where p(r) = and s z (T) = - p(T) taxsh(KW(T)) 



ta,nh(KW(T)) 
KW{T) 



V(r) 



III. MONTE CARLO SIMULATIONS 

These various thermodynamical quantities are determined by Monte Carlo sampling, with respect to the Boltzmann 
wcig ht e" /3W <= tf /Z(T,0). The simulations were first performed using the metropolis method |l7| ]. This method has 
been proved to be very useful at high temperature, far from the phase transition. However, it suffers from a severe 
slowing down near the phase transition, and therefore becomes rather inefficient at low temperature, since, for an 
isotropic two-dimensional system, the critical temperature is Tc = K |is{ ] . A cluster flipping method developed by 
Wolff [fill drastically reduces the slowing down, and the Wolff algorithm was used for high values of K (i.e. low values 



x The numerical values of (P) n , „ and {R) n te as a function of the temperature for the four lattices can be provided 



on request. 
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of T) . To overcome the finite size effects in the critical region, we increased the size of the systems as T decreased. 
The number of cells is 2 14 for K < 2, and increases up to 2 16 for K = 5. Periodic boundary conditions were imposed 
to the system, and a random spin configuration was taken as an initial spin configuration. The number of Metropolis 
steps necessary to reach the termal equilibrium was found to be between 10 2 and 10 4 lattice sweeps depending on the 
temperature. The averaging of the various observables was stopped when A^/x < 0.01. The relative uncertainty on 
the energy and the specific heat was then better than 10~ 2 . These calculations were performed on a Cray J916. 
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FIG. 2. \T (in cm 3 K mol -1 units) as a function of temperature (in units of J). The local spins are assumed to be 
Sc = S'(Mn) = 5/2 and Sq = S(Cu) — 1/2, and the local Zeeman factors g\ = g% = 2. The labels refer to the lattices of figure 
1. 



A. The Magnetic Susceptibility 

Figure 2 shows the behaviour of yT (in cm 3 K mol -1 units) as a function of the temperature (in units of J) for the 
four lattice configurations. The four curves present the same general trend characterised by the following features: 

• a constant value at high temperature corresponding to the paramagnetic limit; 

• a shallow minimum, characteristic of a ferrimagnetic system with antiferromagnetic couplings. It is due to the 
local ordering appearing as the temperature is lowered, and which causes a decrease of the local magnetisation; 

• a rapid increase at low temperature due to the critical divergence at T = 0. 

It is interesting to determine the extent to which all these results can be described by the same universal trend, simply 
corrected by the lattice effects. Actually, at high temperature, the susceptibility is given by the Curie constant which 
depends in a simple way on the spin arrangement: 



35 9 3 9 



(4) 



where uq and ric are the numbers of low spin (quantum) and high spin (classical) ions per molecule, respectively 
(these numbers are given in table I for each lattice). In equation (||) we have taken = 1/2, S 1 -^ = S(Mn) = 5/2. 
For the curves of figure 2 we set gi = gi = 2. 
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TABLE I. Coefficients of the high temperature expansion for c v for the four lattices. 
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The temperature where the minimum occurs is measured from the Monte-Carlo data. We present the results in table 
I (column 4) for each lattice. We have computed this quantity in the mean field theory, and we found the surprinsingly 
simple result : T£(MF) — ^JS 2 n, where n is the number of classical spins in a unit cell (or equivalently the number 
of classical neighbours of a quantum spin). These values are presented in the 5th column of table I. Although the 
mean field and Monte-Carlo results are quantitatively different, one can see from table I that the relative position of 
the minima with respect to the lattice is the same for the two results. This is an indication that the parameter n is 
probably relevant to this quantity. 

At low temperature, the critical divergence is described by the non linear sigma model. The mapping between the 
discrete model on a specific lattice and the universal field theory can be established as follows. First, take the low 
temperature spin wave limit of the lattice model. 



Heff(J) — > Hsw(J* 



<ij> 



Then take the long wave length limit in order to go to the continuum model 

n sw (J*) -^Ha(J) = ij J d 2 r£9„n(r)-^n(r) 



(5) 



The first step depends on the effective interaction between the classical spins whereas the second is related to the 
geometry of the classical spin lattice. The expression of J is given in table I for each lattice model. As a consequence, we 
expect a universal behaviour of the low temperature regime of all the lattices, provided the temperature is renormalised 
in such a way that T = T / J . We do not observe this behaviour quantitatively on our Monte-Carlo data, since our 
lowest temperatures do not lie within the universal critical regime^. However, the hierarchy of the J values gives 
qualitatively well the relative positions of the critical increase of \F for each lattice. 




FIG. 3. The specific heat c„ = 



N c k E 



versus K = 



for the four lattices. 



B. The Specific Heat 

The specific heat is plotted as a function of K in figure 3 for the four types of lattices. In all cases, it presents 
a well pronounced maximum. Beside this common trend, the dependence on the lattice geometry of the details of 
these curves can be understood from simple arguments. The T = (large K) limiting value can be obtained from the 
magnon contribution to the hamiltonian of eq. (^J) which gives an energy per classical spin 



2 we noticed in ref [jl2| that this behaviour is observed for lattice (c) for K > 2.5 which corresponds to T < 0.59J 
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E ^ E + k B T 



Therefore, by normalising the heat capacity to the number of classical spins, we get for T — 0, cy = Cv /Neks = 1 
as can be seen on figure 3 from the result of the simulation. For large T (small K) 1 the behaviour of the specific heat 
can be inferred from the high temperature expansion. The first few terms of this expansion can easily be derived and 
we get 

c v = a 2 K 2 + a A K i + a 6 K 6 H 

where ai turns out to be the number of unit cells connected to a single (classical) site. This number together with 
the coefficient a A and ag are given in table II for each lattice The result, presented in figure 4 in comparison with the 
Monte-Carlo data, shows that up to K ~ 0.4, the system is driven by its high temperature behaviour. The behaviour 
at intermediate and low temperature results on the superposition of two contributions, 
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FIG. 4. Same as figure 3 but in the high temperature regime. The data points correspond to the Monte-Carlo results and 
the lines to the low order high temperature expansion. 

similarly to the case of lattice (c) already analysed in ref.( Jl^]) and where this effect is clearly visible: 

• a first bump for K w 1 resulting from the local ordering of the quantum spin with respect to their randomly 
distributed classical neighbours. By neglecting the correlations between the classical spins, we can estimate this 
contribution 



c$ ~ N Q k B K 2 1 n 



dfi< 



iF 2 (r) 



ier 4tt cosif (KW(T)) 



(6) 



• a second bump at higher value of K, connected to the critical behaviour of the system and which is expected 
to be described by the low temperature limit of the model (eq.(0)). 
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TABLE II. Coefficients of the high temperature expansion for c v for the four lattices. 
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FIG. 5. The specific heat substracted from the quantum spin contribution as a function of the lattice renormalised coupling 
K = jK, (see table I). 

By substracting the specific quantum contribution given by eq.(^) and renormalising the temperature dependence of 
the residual specific heat by the factor defined in the preceding section, we should obtain a universal curve. This is 
what we can observe in figure 5, where the dependence of the lattice has been washed out by comparison with figure 
3. 



IV. CONCLUSION 



In this paper we have determined the effect of the topology of the spin lattice on the thermodynamical properties 
of two-dimensional systems with alternating quantum-classical spins, modeling a wide family of magnetic molecular 
compounds. Since the quantum spin dependence can be traced out, we can used the very powerful classical Monte- 
Carlo techniques to analyse these systems. 

We found that the temperature dependence of both the magnetic susceptibility and the specific heat can be described 
on the ground of general behaviours on which the lattice influence is explicited. This analysis allows one to predict 
in principle what would be the specific heat and the magnetic susceptibility for other compounds with different spin 
geometry. Alternatively, it can be used to determine the interaction parameters together with the Zeeman factors out 
of experimental data pj . 

In most cases |$J, one observes experimentally at low temperature (typically 10 to 15K) a transition towards a 
ferromagnetic ordered phase which cannot be accomodated within our isotropic interaction uM. Therefore, the 
present analysis is only valid in the paramagnetic phase. At low temperature, spin anisotropy and/or spatial 
anisotropy must be taken into account to explain this ferromagnetic phase transition. 
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